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ABSTRACT 

Many asteroids that make close encounters with terrestrial planets are in a binary configuration. 
Here we calculate the relevant encounter timescales and investigate the effects of encounters on a 
binary's mutual orbit. We use a combination of analytical and numerical approaches with a wide 
range of initial conditions. Our test cases include generic binaries with close, moderate, and wide 
separations, as well as seven well-characterized near-Earth binaries. We find that close approaches 
(<10 Earth radii) occur for almost all binaries on 1-10 million year timescales. At such distances, our 
results suggest substantial modifications to a binary's semi-major axis, eccentricity, and inclination, 
which we quantify. Encounters within 30 Earth radii typically occur on sub-million year timescales 
and significantly affect the wider binaries. Important processes in the lives of near-Earth binaries, 
such as tidal and radiative evolution, can be altered or stopped by planetary encounters. 

Subject headings: minor planets, asteroids: general - minor planets, asteroids: individual (2000 DP107, 
1999 KW4, 2002 CE26, 2004 DC, 2003 YTl, Didymos, 1991 VH) 



1. INTRODUCTION 

Near-Earth asteroids (NEAs) live in a dynamic en- 
vironment that includes gravitational encounters with 
planets and the Sun as well as other non-gravitational 
perturbations, with average dynamical lifetimes on the 
order of a few million years (jBottke et al .|l2002f). They 
are replenished from source regions in the main belt of as- 
teroids, which include strong resonances wi th Jupiter and 
Saturn . R adar and ligh t curv e studies bv iMargot et al.l 
(|2002f ) and lPravec et al.l ()2006l ) have determined that ap- 
proximately 15% of these short-lived objects larger than 
200 m in diameter are in a binary configuration, with a 
primary and secondary orbiting their common center of 
mass. The mutual orbits of these binary asteroids are 
affected by tides, by external perturbations such as the 
radia tive binary YORP effect ()Cuk fc BurnsI [20051 : [CuS 
I2007'). and by close encounters with terrestrial planets. 
In this work, we investigate the effect of close planetary 
encounters on the mutual binary orbits of NEA systems 
by studying changes in semi-major axis, eccentricity, and 
inclination. We also calculate the frequency of such ter- 
restrial planet fiybys. 

Early work on the subject of plan etary encounter s with 
NEA binaries was accomplishe d bv iFarinellal (1992) and 
iFarinella fc Chauvineaul (|l993D . They employed analyt- 
ical estimates and Monte Carlo techniques to estimate 
the change in orbital energy and a ngular momentum due 
to an encounter. Recent work by iFang et aD ()2011| ) ex- 
amined the effect of close planetary encounters on the 
mutual orbits of NEA triples 2001 SN263 and 1994 CC. 
They found that scattering events by terrestrial planets 
can excite the eccentricities as well as mutual inclina- 
tions of the satellites' orbits to currently observed val- 
ues on million-year timescales. During close approaches 
with Earth, the outer satellites in both triple systems 
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can have their orbital eccentricities excited to values of 
at least 0.2 as far away as encounter distances of ^^40 
Earth radii for 2001 SN263 and -50 Earth radii for 1994 
CC. Since the orbital effects of encounters on triples have 
been previously examined, this present study focuses 
only on binary systems. Other recent studies have in- 
voked planetary enc ounters to exp lain the rotational dy- 
nami cs of asteroids (jScheeres et al. 2004; iSharma et al.l 
l2006f ). the tidal disru ption of r ubble pile asteroids to 
form binaries (Wal sh fc Richardso n 2006), and t he re- 
surfa cing of NEAs IfNesvorny et al.l 2010: Binzcl et al.l 



20101). Prior stud i es inc lude iBottke fc Mclosh (1996a 
Asphaug fc Bend (|1996f ). and iRichard son et al.l (|199i 



In this work, we perform N-body simulations and em- 
ploy analytical expressions to evaluate a binary's orbital 
changes due to a close planetary encounter with Earth. 
Initial conditions include 3 generic cases (a "close binary" 
with a separation of 4 primary radii, a "moderately- 
separated binary" with a separation of 8 primary radii, 
and a "wide binary" with a separation of 16 primary 
radii) as well as parameters drawn from known, well- 
characterized NEA binaries. By "well-characterized" we 
mean a binary for which the system mass, semi-major 
axis, eccentricity, and approximate component sizes are 
known - in practice this corresponds to a subset of the 
radar-observed binaries. Table [T] shows a compilation 
of these binaries and their parameters such as primary 
size, primary mass, and primary-secondary separation. 
In Section [2l we perform single flyby simulations and 
use analytical equations to determine the orbital effects 
of a close planetary encounter, by thoroughly examin- 
ing a variety of encounter geometries, distances, and ve- 
locities. We study the changes in the semi-major axis, 
eccentricity, and inclination of the binary orbit. In Sec- 
tion [3l we present long-term simulations with test par- 
ticles and planets to calculate encounter timescales for 
observed NEA binaries, whose encounter frequencies are 
strongly dependent on their individual evolutionary his- 
tories. We discuss and summarize this study in Section 
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Table 1 

Well-Characterized Near-Earth Binaries 



System 


Rp (km) 


Mp 


(kg) 


a (km) 


a/Rp 


ao (AU) 


60 


10 (deg) 


(185851) 2000 DPIOT^ 


0.40 


4.38 


X 


10" 


2.62 


6.6 


1.37 


0.38 


8.67 


(66391) 1999 KW4^ 


0.66 


2.35 


X 


1012 


2.55 


3.9 


0.64 


0.69 


38.89 


(276049) 2002 CE26^ 


1.75 


2.17 


X 


1013 


4.87 


2.8 


2.23 


0.56 


47.31 


2004 DC^ 


0.17 


3.57 


X 


lOlO 


0.75 


4.4 


1.63 


0.40 


19.45 


(164121) 2003 YT1_^ 


0.55 


1.89 


X 


1012 


3.93 


7.1 


1.11 


0.29 


44.06 


(65803) Didymos^ 


0.40 


5.24 


X 


10" 


1.18 


3.0 


1.64 


0.38 


3.41 


(35107) 1991 VH^ 


0.60 


1.40 


X 


1012 


3.26 


5.4 


1.14 


0.14 


13.91 



Well-characterized near-Earth binaries and their parameters are compiled in this table, including 
the radius of the primary Rp, mass of the primary Mp, semi- major axis a (in units of km), and 
semi-major axis divided by the primary radius a/Rp. Rough uncertainties in binary parameters 
are ~20% in sizes, ~10% in masses, and <5— 10% in semi-major axes. Heliocentric orbital data 
are given for the semi- major axis aQ, eccentricity bq, and inclination 10 with respect to the 
ecliptic. 

l^Margot et al.' (200T) 

lOstro ct al. (2006) 

tShcDard ct al. (2006) 
Taylor et al. (2008) 

tNolan et al. (2004) 

IBenner et al. (2010) 
Margot et al.i (120081) : IPravec et al.1 l l200ll 

2. EFFECT OF A SINGLE PLANETARY ENCOUNTER 
2.1. Methods 

We perform N-body numerical integrations to inves- 
tigate the effects of a close approach to an Earth-mass 
planet on a binary's mutual orbit. We use a Bulirsch- 
Stoer algorithm fr om an N-body numerical integration 
package, Mercury ()Chamberslll999f ). Our simulations in- 
corporate 3 massive bodies (all assumed to be spheri- 
cal) , which include an asteroid binary and an Earth- mass 
perturber on a hyperbolic trajectory at various close en- 
counter distances. We explore a wide range of initial con- 
ditions in encounter distance and velocity, which we de- 
scribe in detail below. For each pair of encounter distance 
and velocity, systematic simulations are performed with 
nearly 7,000 permutations of the following mutual orbital 
parameters: inclination, longitude of the ascending node, 
and mean anomaly. Our test cases include all of the NEA 
binaries shown in Table [T] as well as a few generic cases 
with a typical rubble pile density of 2 g cm~^ and a pri- 
mary radius of 0.5 km: a "close binary" with a separation 
of 4 primary radii, a "moderately-separated binary" with 
a separation of 8 primary radii, and a "wide binary" with 
a separation of 16 primary radii. 

The hyperbolic encounter velocity Voo, defined as the 
relative speed between the binary and Earth at infinity, 
can be described by a probability distribution of relative 
velocities as shown in Figure [1] As a result, our sim- 
ulations appropriately cover a range of Voo from 8—24 
km s~^ with increments of 4 km s~^. We also systemati- 
cally repeat the ensemble of simulations for various close 
encounter distances whose range is dependent on the bi- 
nary's semi-major axis, and with typical increments of 
1-2 

In this paper, the "encounter distance" q is defined 
as the distance of closest approach and the "impact pa- 
rameter" b is the hypothetical encounter distance that 
would result in the absence of gravitational focusing. 
The impact parameter can be described as the perpen- 
dicular distance between the non-focused path of the 
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10 20 30 40 

Encounter Velocity (km/s) 

Figure 1. This probability distribution of expected en- 
counter velocities Woo of NEAs with Earth is obtained 
from numerical simulations of asteroid migration from 
source regions (Mars-crossing region, z/g resonance, and 
3:1 mean-motion resonance with Jupiter) in or adjacent 
to the main belt (Bill Bottke and Kevin Walsh, personal 
communication, 2011). 

perturber and the binary. The usual definition relat- 
ing impact parameter b and encounter distance q is 
62 ^ q2(^i (2G'Me)/((7u^)), where G is the gravita- 
tional constant, is the mass of Earth, and v^c is the 
encounter velocity at infinity. For most cases we consider 
here, the difference between b and q is less than 1 Earth 
radius. 

All of the simulations begin with asteroid binaries in 
circular orbits. After a scattering encounter, possible 
outcomes include (a) an intact, stable binary system, 
or instability marked by (b) collision between any two 
bodies or (c) ejection of the secondary from the binary 
system. Encounter results are recorded for changes in 
semi-major axis, eccentricity, and inclination in stable 
encounters only. Stable systems are defined as binaries 
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with no collisions nor ejections of the secondary. 

The output from numerical simulations is compared 
with approximate analytical results. We consider close 
planetary encounters with NEA binaries to be impulsive 
events, defined as swift encounters where the planet's 
interaction with the binary is much shorter than the or- 
bital period of the binary. The encounter delivers an 
impulse representing a shift in velocity of the binary's 
orbit, and this impulse approximation is valid in the do- 
main where q/voa *S where P is the binary's mu- 
tual orbital period. For such impulsive encounters, the 
change in a binary's orbital elements h as been analyti- 
cally a pproximated by previous studies. iHeggie fc Rasiol 
(1 19961) derived the change in a binary's eccentricity, 
iCoUins fc Saril ()2008l ) derived the chan ge in both ec- 
centri city a nd inclination, and [Pa rinell a fc Chauvineaul 
(|1993l ) and iChauvineau fc Farinellal (11995^ derived the 
change in energy (which we can relate to a change in 
semi-major axis). The typical change in a binary's or- 
bital elements after averaging over all encounter geome- 
tries (which include stable and unstable encounters) can 
be estimated as 



Aa « 1.48 



Ae 



Ai « 0.75 



G Mffia5/2 




(1) 



(2) 



(3) 



where G is the gravitational constant, is the mass 
of Earth, Vao is the encounter velocity, and q is the en- 
counter distance. The binary's semi-major axis is repre- 
sented by a, its eccentricity is e, its inclination is i, and its 
system mass is M . The constants in front of Equations 
[2] and [3] represent the averagi ng of encounter geome tries, 
and have been calculated by iCoUins fc Saril ()2008[ ). In 
their paper, we note there is a factor of 2 missing in the 
denominator of Equation A3. By analogy, we obtain the 
constant in front of Equation [T] by scaling the analytical 
curve to our numerical simulation results. 

2.2. Results 

The change in semi-major axis due to a planetary flyby 
is shown for three generic cases presented in Figure [2l 
Output from numerical simulations are shown as points 
and represent the mean values of the change (after tak- 
ing the absolute value) in semi-major axis resulting from 
stable encounters. The curves represent analytical esti- 
mates, which provide a reasonable match to our numeri- 
cal results. We show the averages changes in semi-major 
axis for Voo = 12, 16, 20, and 24 km s~^. The results 
using a Voo of 8 km s~^ are not shown because the re- 
sulting spread in the post-encounter semi-major axis is 
very broad (e.g. at an encounter distance of 10 Earth 
radii, the standard deviation is ~2.1 km) and not ac- 
curately portrayed by a single value. We wish to show 
orbital effects of only stable encounters; the range of close 
encounter distances for each binary type shown in this 
figure is determined by the distances at which all encoun- 
ters resulted in stable binaries. Since the 3 types of bi- 
naries have different separations, they will have different 



ranges of encounter distances at which stable encounters 
occur. These plots show that larger binary separations 
produce greater changes in semi-major axis at any given 
encounter distance because wider binaries are less tightly 
bound and thus more susceptible to passing perturbers. 
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Figure 2. The change in a binary's semi-major axis is 
shown as a function of encounter distance and velocity for 
3 types of binaries: close binaries, moderately-separated 
binaries, and wide binaries. Results from numerical sim- 
ulations are shown as dots and analytical calculations are 
depicted by solid lines. 

The change in eccentricity is shown in Figure [Sj with 
similar encounter distance and velocity ranges as shown 
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Figure 3. The change in a binary's eccentricity is shown 
as a function of encounter distance and velocity for 3 
types of binaries: close binaries, moderately-separated 
binaries, and wide binaries. Results from numerical sim- 
ulations are shown as dots and analytical calculations are 
depicted by solid lines. 

in Figure [21 Results from numerical simulations are 
given as the mean value of the eccentricity increase of 
all stable encounters. Analytical and numerical results 
are in agreement, especially for greater encounter veloc- 
ities. For our test NEA binaries (Table [T]) and using a 
typical Voo of 12 km s~^, we also show the close encounter 
distances at which their eccentricities can be excited to 
a value of at least 0.2 (Figure HJ and how the distance 
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Figure 4. Using Orbital and physical parameters from 
the sample of NEA binaries given in Table [U we show 
the fraction of stable systems (no ejections nor collisions) 
with excited eccentricities (e > 0.2) as a function of en- 
counter distance for v^o = 12 km s~^. 
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Figure 5. The eccentricity excitation distance (defined as 
the close encounter distance where 50% of stable NEA 
binary systems showed e > 0.2) is given as a function of 
the secondary's semi-major axis in units of primary radii. 
The curve shows the result obtained from the analytical 
expression. This case is for Uoo of 12 km s^^. 

at which eccentricity excitation occurs is a function of 
binary separation (Figure [5]) . The eccentricity excitation 
distance is defined in Figure [5] as the encounter distance 
where 50% of stable encounters resulted in binaries with 
e > 0.2. Since it is evident from this figure that there 
is a relationship between the critical encounter distance 
q and the binary's semi-major axis a (when expressed in 
units of primary radii), we solve for its analytical rela- 
tionship by rearranging Equation [2l and find 



g 



1 



\A7rp) 



Rr 



3/4 



(4) 



where Rp is the primary's mass, i?^ and are the 
radius and mass of Earth, G is the gravitational constant, 
Ae is the change in eccentricity due to the encounter, Voo 
is the encounter velocity, and p is the binary's density. 
As stated earlier, the 1.89 factor is due to averaging 
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Figure 6. These plots show the characteristic behav- 
ior of close NEA binaries (i.e. 1999 KW4), moderately- 
separated NEA binaries (i.e. 2000 DP107), and wide 
NEA binaries (represented by a hypothetical binary with 
a separation of 16 primary radii) due to varying en- 
counter velocities {voo — 8—24 km s^^ with increments 
of 4 km s~^). The fraction of stable systems (no ejec- 
tions nor collisions) with excited eccentricities (e > 0.2) 
is shown as a function of encounter distance. 



over all angles of encounter geometries ([Collins fc Saril 
I2008f ). This equation is valid for initially circular bina- 
ries. The same q oc a^^^ relationship holds for all bi- 
naries in our sample because for a uniform set of en- 
counter parameters (encounter velocity and perturber's 



radius and mass), encounter strength (strong enough to 
excite Ae = 0.2), and binary density (assuming NEA bi- 
naries have similar rubble pile densities), the pref actor in 
front of the a^^^ term is constant. In convenient units, 



R(Tl 



3.21 



0.2 12 km s-i 



1/2 



2 g cm 
P 



1/4 



Rp 
(5) 

We overplot this analytical relationship in Figure [5] 
and find a good match to the numerical results. The 
agreement between numerical and analytical estimates 
for 2003 YTl is not as good due to its higher density 
estimate (its nominal density is p 2.7 g cm^^), since 
we assumed p — 2 g cm^'^ in the calculation of Equation 

El 

In Figure [S] we show the eccentricity excitation be- 
havior due to VoD = 8, 12, 16, 20, and 24 km s~^ 
for specific observed NEA binaries: a close binary (de- 
picted by 1999 KW4), a moderately-separated binary 
(2000 DP107), and a wide binary (represented by a hy- 
pothetical binary with a 16 primary radii separation). 
As the value of Uoo decreases, a greater fraction of sta- 
ble systems have excited eccentricities due to a planetary 
flyby at a given encounter distance. In addition, as the 
value of Voo decreases, the final eccentricities of binaries 
in our simulations increased. These effects occur because 
a slower-passing perturber will have a longer encounter 
duration with the binary and therefore cause a stronger 
perturbation on the binary's mutual orbit. 

The change in a binary's orbital inclination is shown 
in Figure [71 with similar encounter distances and veloci- 
ties as shown in Figures [2] and [3l Results from numerical 
simulations are given as the mean value of the inclination 
change of all stable encounters. Comparison between 
output from numerical simulations and analytical esti- 
mates provide a decent match. Possible reasons for the 
less-than-exact agreement between numerical and ana- 
lytical results include our choice of using the mean value 
to represent the average change in inclination of stable 
encounters from simulations as well as the practical limit 
on the number of encounter geometries we were able to 
perform in our simulations. The analytical expressions 
are also only valid in the limit of impulsive encounters, 
suggesting that analytical results should provide better 
matches to simulations with faster encounter velocities, 
which is consistent with our results. 

These results have implications for binaries observed 
with excited dynamical states. For instance, previous 
studies of 1999 KW4 (heliocentric a is 0.64 AU, e is 0.69, 
i is 38.89 degrees; current solar close approach distance 
is 0.2 AU) suggest that the binary's close approaches 
to the Sun can excite its rotational an d orbital states 
([Ostro et all 120061 : [Scheeres et all 120061 ). Their numeri- 
cal simulations show that the binary's mutual orbit pole 
can be modified by more than 0.5 degrees per pericen- 
ter passage for a pericenter distance of 0.2 AU and by 
more than 1 degree for a pericenter distance of 0.12 AU. 
From our simulations, we find that planetary encounters 
with an Earth-mass body can perturb the binary orbit 
pole by a similar amount. Therefore, even NEAs that do 
not come as close to the Sun as 1999 KW4 may exhibit 
excited rotational and orbital dynamics. For a generic bi- 
nary with a semi-major axis of 4 primary radii (same as 
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Figure 7. The change in a binary's inchnation is shown 
as a function of encounter distance and velocity for 3 
types of binaries: close binaries, moderately-separated 
binaries, and wide binaries. Results from numerical sim- 
ulations are shown as dots and analytical calculations are 
depicted by solid lines. 

1999 KW4's separation), our simulations using Voo from 
8—24 km s~^ show that a 0.5 degree shift in inclination 
can occur due to a planetary encounter at distances of 
18—28 i?0 and a 1 degree shift can occur at 13—20 R^. 

We briefly discuss instability trends seen in simulations 
with a typical v^o of 12 km for observed NEA bina- 
ries (Table[T]). For nearly all encounter distances sampled 
in simulations, ejections dominated over collisions in un- 



stable encounters. If planetary encounters close enough 
to disrupt a binary have occurred, an ejected secondary 
is more likely to occur than a collision between the pri- 
mary and secondary. Therefore one would expect that 
planetary encounters can lead to the formation of aster- 
oid pairs in the near-Earth population, in addition to 
the formation mechanisms responsible for pair produc- 
tion in the main belt. The different routes towards the 
creation of asteroid pairs are d iscussed in a companion 
paper bv lFang fc Margotl ([20T1 . Close binaries such as 
Didymos, 2002 CE26, 1999 KW4, and 2004 DC have sim- 
ilarly small component separations (2—5 Rp) relative to 
their primary's radius Rp, and exhibit comparable dis- 
ruption statistics. Instabilities, including both ejections 
and collisions, occurred at encounter distances of ~5— 7 
i?© and less. More moderately-spaced binaries such as 
2000 DP107, 2003 YTl, and 'l 991 VH have compara- 
ble component separations (5—9 Rp). For these bina- 
ries, simulation results showed disruption occurring at 
encounter distances of ~7— 9 i?© and less. Binaries with 
wide separations will have mutual orbits that are more 
easily perturbed than tighter systems. 

Given the significant orbital perturbations in binary 
systems that can be caused by planetary flybys, the next 
section examines the frequency of such encounters. 

3. FREQUENCY OF PLANETARY ENCOUNTERS 
3.1. Methods 

After investigating the orbital effects of a single plan- 
etary flyby on a binary in Section [21 we next determine 
the frequency of such close planetary encounters with 
Earth for well-characterized NEA binaries (Table[T]) . The 
timescale of close encounters is very dependent on the 
specific NEA's past evolutionary path to its current helio- 
centric orbit (see Table[l]for current heliocentric parame- 
ters), and different NEAs can have very different close en- 
counter histories. We do not backward-integrate the cur- 
rent heliocentric orbits because the orbits become chaotic 
over time and it is difficult to reconstruct their past 
history, eve n statistically . Con sequently, we follow the 
methods of IBottke et al l (j2002l) by integrating test par- 
ticles from the three strongest source regions {ve secular 
resonance with Saturn, 3:1 mean-motion resonance with 
Jupiter, and Mars-crossing regions) in or adjacent to the 
main belt and by tracking the test particles as they mi- 
grate into near-Earth space. These source regions are the 
main dynamical pathways to NEA orbits. Integrations 
are performed using a hy brid symple ctic/Bulirsch-Stoer 
algorithm from Mercury (jChambersI [1999) . and include 
the Sun and all 8 planets in addition to 3000 test parti- 
cles per source (9000 test particles total). In the simula- 
tions, we ignore the effect of Yarkovsky. The Yarkovsky 
effect is an important perturbation that moves asteroids 
into main source regions; however, once asteroids have 
reached the source regions, the effect of Yarkovsky is neg- 
ligible compared to the effect of strong resonances and 
planetary encounters. 

After test particles have migrated into near-Earth 
space, we search for particles whose orbital elements 
(semi-major axis a, eccentricity e, and inclination i) 
closely match the current orbital elements of the actual 
NEA binaries in our sample (Table [ij. For each bi- 
nary, we search for ten test particles closest in orbital 
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Table 2 

Encounter Timcscales (Myr) with Earth 



Name 


Impact Parameter: 

2R(s 4Re 6i?e 


8R(B 


lOiJe 


12R© 


I4i?e 


16iJe 


18R(B 


20R(B 


22iJ® 


24R(s 


26i?© 


28Re 


30i?e 


2000 DP107 


60.73 


15.18 


6.75 


3.80 


2.43 


1.69 


1.24 


0.95 


0.75 


0.61 


0.50 


0.42 


0.36 


0.31 


0.27 


1999 KW4 


50.72 


12.68 


5.64 


3.17 


2.03 


1.41 


1.04 


0.79 


0.63 


0.51 


0.42 


0.35 


0.30 


0.26 


0.23 


2002 CE26 


391.40 


97.85 


43.49 


24.46 


15.66 


10.87 


7.99 


6.12 


4.83 


3.91 


3.23 


2.72 


2.32 


2.00 


1.74 


2004 DC 


109.90 


27.47 


12.21 


6.87 


4.40 


3.05 


2.24 


1.72 


1.36 


1.10 


0.91 


0.76 


0.65 


0.56 


0.49 


2003 YTl 


131.34 


32.84 


14.59 


8.21 


5.25 


3.65 


2.68 


2.05 


1.62 


1.31 


1.09 


0.91 


0.78 


0.67 


0.58 


Didymos 


46.30 


11.58 


5.14 


2.89 


1.85 


1.29 


0.94 


0.72 


0.57 


0.46 


0.38 


0.32 


0.27 


0.24 


0.21 


1991 VH 


72.86 


18.21 


8.10 


4.55 


2.91 


2.02 


1.49 


1.14 


0.90 


0.73 


0.60 


0.51 


0.43 


0.37 


0.32 



For each NEA binary, we show the encounter timescale (in millions of years) with Earth for various impact parameters. 



element space to the binary's current heliocentric or- 
bital elements using this scaling : 5a/ a ^ (5e ^ (5sini 
(|Vokrouhlickv fc Nesvorn?! |2008| ). Binaries 2003 YTl 
and 1999 KW4 have either high inclinations and/or low 
semi-major axes, and these binaries proved to be difficult 
when searching for matches with test parti cles originat- 
ing fro m the 3 main source regions given by lBottke et alJ 
(j2002f ). In our sample of binaries, the median da is 
~0.014 AU. The matching test particles are further inte- 
grated until they become dynamically unstable through 
ejection from the Solar System or collision with a planet 
or the Sun. For each NEA binary's matching test parti- 
cles, we record the orbital history with a n Opik-type code 



similar to the proce dure developed by IWetherilll ()1967l ) 
and implemented by iFarinella fc Davia (119921 ) . which al- 
lows us to calculate the encounter velocity and intrinsic 
coUisional probability with Earth at each timestep. We 
calculate the intrinsic collisional probability with Earth 
every 10,000 years from the time of injection in the source 
region up to the most recent epoch at which the test 
particle's orbital elements matched those of the actual 
NEA binary. As a result, we obtain an intrinsic col- 
lisional probability value for each 10,000-year timestep 
in the history of each matching test particle. Based on 
this ensemble of values we associate an average intrinsic 
collisional probability to each test particle by averaging 
all non-zero values over time. Finally, we associate an 
intrinsic collisional probability to each NEA binary by 
averaging over all of its matching test particles. 

The intrinsic collisional probability (in units of km^^ 
yr^^) can be multiplied by the square of an impact pa- 
rameter to obtain an encounter probability (in units of 
yr~^), whose inverse yields an approximate encounter 
timescale. These timescales represent the average time 
between encounters with Earth since injection into main 
belt source regions and their subsequent migration into 
near-Earth space. 

We also calculate the close encounter timescale 
using an analytical, order- of-magnitude approach. 
iChauvineau fc Farinellal ()1995D give the approximate 
timescale t 

(6) 



t 



between close encounters up to an impact parame- 
ter b given an asteroid's lifetime Tcoh against colli- 
sion with Earth of radius i^m. For a given asteroid 
size, iStuart fc Binzell (|2004[ ) provide the total number 



of NEAs with at least that size as well as the average 
interval between Earth impacts due to all asteroids of 
that size taken collectively. Thus, with these numbers 
we can calculate Tcoh and the encounter timescale for an 
individual asteroid. 

3.2. Results 

For each binary in our sample. Table gives encounter 
timescales with Earth from numerical simulations for im- 
pact parameters ranging from 2Rq, to 3Oi?0 . In this im- 
pact parameter range, encounter timescales can range 
from '--^10^ up to '^10^ years for the NEA binaries con- 
sidered here. Encounter probabilities can be obtained by 
taking the inverse of these timescales, and then suitably 
scaled (the probability varies as the square of the impact 
parameter) to a desired impact parameter. Recall that 
due to gravitational focusing, the impact parameter b is 
not the same as the closest encounter distance q, and 
are related in this manner: 5^ = q'^{l + (2GMQ)/{qvl^)), 
where G is the gravitational constant, is the mass of 
Earth, and Uqo is the encounter velocity at infinity. 

Timescale estimates from the analytical approach 
agree with results obtained from numerical integrations 
within an order of magnitude for all NEA binaries in our 
sample. Differences in timescales can be attributed to 
the fact that the analytical timescales do not take into 
account a specific NEA's orbital trajectory as it migrates 
to near-Earth space, which affects its history of planetary 
encounters. 

We repeat this analysis to determine encounter fre- 
quency with other terrestrial planets, namely Mercury, 
Venus, and Mars. For a given orbital effect such as in- 
creasing a binary's eccentricity to a given value, we find 
that encounters with Venus are just as important as en- 
counters with Earth. For all binaries in our sample, the 
encounter probabilities with Venus and Earth are compa- 
rable; these encounter probabilities for achieving a given 
change in the mutual orbit are a function of encounter 
velocity, planet mass, and binary separation. 

4. DISCUSSION AND CONCLUSION 

Encounter timescales presented in Table [2] show that 
some NEA binaries in our sample (Table [1]) can typically 
encounter Earth frequently enough at close-enough dis- 
tances to excite the binary's orbital elements, as shown 
by single-fiyby simulations in Section [21 We have pre- 
sented the encounter timescales in this study, but the 
actual effects of encounters are also dependent on typi- 
cal encounter velocities, which we plotted in Figure[T]for 



7 



a generic NEA. Consideration of these factors as well as 
comparison to obser ved eccentricities are d iscussed in a 
companion paper bv lFang fc MargotI pOlTI ). 

We briefly address the possibility of repeat encounter 
passes. Repeat passes, which are most relevant for bi- 
naries with short encounter timescales and large binary 
separations, can increase or decrease the binary's semi- 
major axis and eccentricity depending on its values prior 
to the encounter. The strength of repeat passes depends 
on flyby parameters such as encounter velocity. The net 
effect of repeat passes is still an eccentric orbit, since it 
is rare for an eccentric binary to undergo a planetary 
encounter and end up with near-zero eccentricity. 

Close planetary encounters have important implica- 
tions for an NEA binary's evolution, since flybys can 
disrupt main evolutionary processes in a binary by ex- 
panding or contracting the mutual orbit. For a typical 
NEA binary with a separation of ~4 primary radii, we 
find that the semi-major axis increases on average 60% of 
the time for encounter distances from 2—10 Rq, and en- 
counter velocities from 8—24 km s~^. Since flybys can in- 
crease or decrease the semi-major axis of the binary's or- 
bit, tidal evolution can strengthen (if the semi-major axis 
decreases) or weaken (if the semi-major axis increases) 
since the tidal torque scales as the binary separation to 
the sixth power. When the semi-major axis is modified, 
the mean motion of the secondary changes and an ini- 
tially spin-locked secondary may become asynchronoutQ. 
Asynchronization can also occur by a change in the sec- 
ondary's rotation rate due to the flyby. The loss of spin- 
lock would imply that radiative perturbations depen- 
dent on synchronization such as binary YORP (BYORP) 
would be shut down. Lastly, planetary encounters are 
another mechanism that may form near-Earth asteroid 
pairs. The evolution of binaries under various influences, 
including the planetary model present ed here, is fur- 
ther d iscussed in a companion paper by iFang fc Margoti 
((20T1 . 

To summarize, in this study we have used numerical 
integrations and analytical expressions to investigate the 
effects of planetary encounters on NEA binaries. We 
found the encounter distances at which flybys can in- 
crease the orbital semi-major axis, eccentricity, and in- 
clination for a variety of encounter velocities. There is 
reasonable agreement between results obtained from sim- 
ulations and analytical methods. The possible outcomes, 
including collisions, ejections, and stable encounters, are 
discussed for close binaries, moderately-separated bina- 
ries, and wide binaries. We have also used N-body 
integrations to examine the past evolutionary histories 
of NEAs as they migrated from main belt source re- 
gions into near-Earth space. From these simulations, we 
calculated encounter probabilities and timescales for all 

^ In this paper, binaries with an absence of spin-orbit synchro- 
nism arc called asynchronous binaries. Binaries with a secondary 
spin period synchronized to the mutual orbit period are called syn- 
chronous binaries. Binaries with both primary and secondary spin 
periods synchronized to the mutual orbit period are called dou- 
bly synchronous binaries. Most NEA binaries are synchronous. 
Note that our terminology is different from that of Pravcc & Ifarris 
([2007), who used the term "asynchronous binaries" for binaries 
with spin-orbit synchronization. If generalization to systems with 
more than one satellite is needed, we affix the terms synchronous 
and asynchronous to the satellites being considered. 



NEA binaries in our sample for different impact param- 
eters. These encounter timescales provide rough agree- 
ment with analytical estimates, which do not take into 
account individual NEAs' past orbital histories. Lastly, 
planetary encounters have important implications for 
evolutionary processes such as tidal and BYORP mech- 
anisms. 

We thank Bill Bottke and Ben Collins for useful dis- 
cussions. We are also grateful to the reviewer for helpful 
comments. This work was partially supported by NASA 
Planetary Astronomy grant NNX09AQ68G. 
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